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Abstract 

The absence of the expected GZK cutoff strongly challenges the notion that the highest-energy cosmic rays 

are of distant extragalactic origin. We discuss the possibility that these ultra-high-energy events originate 

in our Galaxy and propose that they may be due to iron nuclei accelerated from young, strongly magnetic 

neutron stars. Newly formed pulsars accelerate ions from their surface through relativistic MHD winds. We 

Q>^ ' find that pulsars whose initial spin periods are shorter than ~ 4{Bs/W^^G) ms, where Bs is the surface 

G^ ■ magnetic field, can accelerate iron ions to greater than lO^'^ey. These ions can pass through the remnant of 

^ ' the supernova explosion that produced the pulsar without suffering significant spallation reactions. Depending 

^ ^ on the structure of the galactic magnetic field, the trajectories of the iron ions from galactic sources can be 

^ ■ consistent with the observed anival directions of the highest energy events. 
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1 Introduction: 

The detection by AGASA (Hayashida 1994; Takeda et al. 1998, 1999), Fly's Eye (Bird et al. 1995, 
1994, 1993), and Haverah Park (Lawrence, Reid & Watson 1971) of cosmic ray events with energies above 
the expected Greisen-Zatsepin-Kuzmin (GZK) cutoff (Greisen 1966; Zatsepin & Kuzmin 1966) has triggered 



(N 

> 

oo 

fr^ , considerable interest (see e.g, Cronin 1999; Bhattacharjee & Sigl 1999 for recent reviews). The cutoff should 

cn ' be present if the ultra-high energy particles are protons from extragalactic sources. Cosmic ray protons of 

^ ■ energies above 5 x 10^^ eV lose their energy to photopion production off the cosmic microwave background 

Q\ [ and cannot originate further than about 50 Mpc away from us. Alternatively, if ultra-high-energy cosmic 

^ ' rays (UHECRs) are protons from sources closer than 50 Mpc, the arrival directions of the events should point 

^ ■ toward their source. The present data shows a mostly isotropic distribution and no sign of the local distribution 
of galaxies or of the galactic disk (Takeda et al. 1998). In sum, the origin of these ultra-high energy particles 



O ' remains a mystery. Here we discuss how the early evolution of young pulsars may be responsible for the yet 

unexplained flux of cosmic rays beyond the GZK cutoff (Blasi, Epstein, & Olinto 1999). 



2 UHECRs from Newborn Pulsars: 



When neutron stars are born during the collapse of supermassive stars, they begin their life with very fast 
r> I rotation (17 ~ 3000s^^) and very large magnetic fields (on the star's surface Bg^W^^ G). Inside the light 

C^ ■ cylinder (i.e, r<Ric where Ric = c/Q,) a magnetosphere corotates with a dipole magnetic field component 

that scales as B{r) = Bs{Rs/fY where the radius of the star Rs ^ 10^ cm. The magnetosphere has a density 
(Goldreich & Julian) nGj{r) ~ i?(r)i7/47rZec , where Z is the charge of the magnetospheric nucleus, e the 
electric charge and c the speed of light. As the distance from the star increases, the dipole field structure cannot 
be maintained and beyond the light cylinder the field is mostly azimuthal. For young rapidly rotating neutron 
stars, the light cylinder is just about ten times the star radius, Ri^. = lO^^ig^^ cm (where 0.^}^ = 17/3000 s^^). 

Iron nuclei stripped by strong electric fields from the surface of the neutron star, fill much of the mag- 
netosphere. From the light cylinder, a relativistic plasma expands as a magnetohydrodynamic (MHD) wind 
(Gallant & Arons 1994). The field strength {B oc r~^) and geometry in this region are such that the plasma 
moves relativistically with Alfven speed close to the speed of light. In the rest frame of the wind, the plasma 
is relatively cold while in the star's rest frame the plasma moves with 7 ^ 10^ — lO^*'. 

The typical energy of the accelerated cosmic rays can be estimated by considering the magnetic energy 
per ion at the light cylinder Ecr = Bf^/Sirncj- At the light cylinder the magnetic field strength is Bic = 



10^° GSiaJ^^fc and ncj = 1.7 x 10^^ cm^^ Bi^Vl\,^/Z which gives 

Eer = 4xl020eVZ26Si3J7L> (1) 

where Z = 26^26 and Bs = 10^'^ G -B13. The spectrum of accelerated UHECRs is determined by the evolu- 
tion of the rotational frequency: As the neutron star ages, the rotation speed decreases due to electromagnetic 
and gravitational radiation. For initial periods of ^1 ms and i?5^10^^ G, the spin down is dominated by 
magnetic dipole radiation given by: IQCi = —B'^RgQ.^/Gc^. 

For a moment of inertia / = 10^^ g cm^, the star's spin down time is tgd = ^/Cl ~ 1.8 x lO^s B^^VLT^^, 
and the spectrum can be found to be (Blasi, Epstein, & Ohnto 1999): 

A,(g),g,$g,g,f^^-'fGeV-. (2) 

Jl dE B1SE20Z2Q 

The particle flux is thus Q = i uqj T^Rfc^ = ^^§-^6 x 10^^ s~^ where ^ < 1 is the efficiency for accel- 
erating particles at the light cylinder. The energy as a function of spin frequency in Eq.Q gives dE/d^ = 
1.7 X 10^3 E/nsk- 

Young neutron stars are most likely distributed in the Galactic disk at typical distances ~ 8 kpc. If they are 
produced at a rate 1/r where r = 100 T2 yr, the flux of particles on Earth is 

F(E)^^^ ('"'" QeV-WV. (3) 

ttR^t T2BisE2oZ2e 

Assuming an isotropic distribution of the arrival directions, the AGASA experiment finds a flux at lO^'^eV of 
F(10^°eV) = 4 X 10~^° GeV~^cm~^s~^. Therefore, the efficiency need only be ^>4 x 10~^. Furthermore, 
the spectrum derived above is very flat, N{E) oc E~^, another property in good agreement with the data. 

Even though the young neutron stars are surrounded by the remnants of the presupernova star, the accel- 
erated particles can easily escape the supernova remnant without significant degradation for a wide range of 
initial magnetic fields and spinning rates. The supernova event that formed the young neutron star also ejected 
the envelope of the original star, making it possible for cosmic rays to escape. However, as the envelope 
expands, the young neutron star spins down and may become unable to emit particles of the necessary en- 
ergy. Thus, an additional requirement for this scenario is that the column density of the envelope becomes 
transparent before the spinning rate of the neutron star decreases significantly. 

To estimate the evolution of the column density of the envelope, consider a supernova that imparts Esn = 
lO^^E^i erg on the stellar envelope of mass M^nv = WMiMq. The envelope then disperses with a velocity 
Ve — {2EsN/Menv) = 3 X 10^ (E^i/Mi) ' cm s~^ . The column density of the envelope surrounding 
the neutron star is S ~ Menv/^T^R^rt where Reff = Ro + Vet. The initial value Rq is characteristic of 
pre-supemova stars and Rq ~ 10^^cmi?i3. We now have 

S ~ 7T = 1.6 X 10 g cm ^ ^. — ^Vt^ ) 

where t is in seconds, and te = Ro/ve = 3 x 10^ Ris\/Mi/E^is. The condition for iron nuclei to traverse 
the supernova envelope without significant losses is that S < Sj ~ 100 g cm~^. At late times compared to 

7 —1/2 

tg, this "transparency" condition gives t > ttr = 1.3 x 10' AdiE^^ s. 

At the same time, the neutron star loses energy at a rate given by dipole radiation whose solution is 



[i + tsBl,nl,] 



where Qi^k is the initial spin period and tg = t/lO^s. The evolution of the maximum energy is thus given by 

Ecr (t) 



4xl020eV ^26i?i3f^L 



[1 + tsBf.nl,] • 

Since the maximum energy decreases as the source evolves, the condition that a source could produce the 
UHECRs is that £„■ exceeds the needed energy when the envelope becomes transparent; i.e., Ecr{ttr) > 
10^'^-E20 eV. This translates into the following condition: 
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^From this equation we obtain upper bounds on the surface magnetic field strength and the star's initial spin 

period. Pi = 2-^/^,,; i.e., Bi^ < 2,IZ2gEI{^ /M1E20 and Pi < ^ttBII^^ Z2&E2Q . For Mi = 2 and E20 = 
E^i = Z2& = 1, this is just i3i3 < 15.4 and Pj<10ms, not very restrictive values for a young neutron star. 
The allowed regions in the Bs-^i plane are shown in Figure 1 for £'20 = 1 and 3. The iron ejected with 
energies ~ lO^'' eV will reach Earth after some diffusion through the Galactic and halo magnetic fields. The 
gyroradius of these UHECRs in the Galactic field of strength Bgai is 
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which is a few times the typical distance to a young neutron star {--^ 8 kpc). Therefore, the ultra-high energy 
iron arriving at the Earth would not point at the source. According to Zirakashvili et al. (1998), a Galactic iron 
source is consistent with the arrival direction distribution observed by AGASA for UHECRs. 

In support of this interpretation, we note that the cosmic ray component at 10^^ eV is nearly isotropic. If 
these cosmic rays are protons of galac- 
tic origin, the isotropic distribution ob- 
served at these energies may be indica- 
tive of the diffusive effect of the Galac- 
tic and halo magnetic fields. The iron 
arrival distribution at 10^*^ eV probes 
similar trajectories to protons at a few 
times 10^^ eV. Depending on the exact 
strength and structure of the magnetic 
field in the Galaxy and the halo of the 
Galaxy, the ultra-high energy iron may 
retain some memory of the source di- 
rection. As the data improves at the 
highest energies a preference for the 
galactic plane should become evident. 
In addition, the relativistic wind may 
also accelerate some lighter nuclei that 
can help constrain the scenario and the 
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Figure 1: Allowed region of parameters for efficient acceleration. 



structure of the Galactic and halo magnetic fields. 

3 Conclusion 

We have discussed the possibility that the ultra-high-energy events observed past the GZK cutoff originate 
in our Galaxy and are due to iron nuclei accelerated from young, strongly magnetic neutron stars. Iron from 



the surface of newborn neutron stars are accelerated to ultra-high energies by a relativistic MHD wind. Pulsars 
whose initial spin periods are shorter than ~ 4(i?5'/10^'^G) ms can accelerate iron ions to greater than lO^'^eF. 
These ions can pass through the remnant of the supernova explosion that produced the pulsar without suffering 
significant spallation reactions. 

This proposal can be tested by future experiments such as the Auger Project. The best test of this proposal 
is a unambiguous composition determination and a correlation of arrival directions for events with energies 
above 10^'^ eV. Both aspects will be testable with future experiments. Finally, a detailed study of the lighter 
element components expected to reach us from this process will also help constrain this scenario. 

Acknowledgments 

The research was partly supported by NSF through grant AST 94-20759 and DOE grant DE-FG0291 
ER40606 at the University of Chicago and partly carried out under the auspices of the U.S. Department of 
Energy with support by IGPP at LANL. 

References 

Bird, D.J. et al., Phys. Rev. Lett. 1993, 71 3401; ApJ 1994, 424, 491; ApJ 1995, 441, 144 

Bhattacharjee, P, & Sigl, G. 1999, Phys. Reps, submitted, astro -ph/98 11011 1 

Blasi, P., Epstein, R.I., & Ohnto, A.V., 1999, to be submitted 

Cronin, J.W., Rev. Mod. Phys. 1999, 71, S165 

Gallant, Y.A., & Arons, J. 1994, ApJ 435, 230 

Greisen, K., Phys. Rev. Lett. 1966, 16, 748 

Hayashida, N., Phys. Rev. Lett. 1994, 73, 3491 

Lawrence, M.A., Reid, R.J.O. & Watson, A. A., J. Phys. G Nucl. Part. Phys. 1991, 17, 733 

Takeda, M. et al, Phys. Rev. Lett. 1998, 81, 1163; ibid, 1999, astro-ph/9902239 submitted to Astrophys. J. 

Zatsepin, G.T, & Kuzmin,V.A, Sov. Phys.-JETP Lett. 1966, 4, 78 

Zirakashvih, V.N., Pochepkin, D.N., Ptuskin, V.S., & Rogovaya, S.I., Astron. Lett. 1998, 24, 139 



